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Preparation of entangled pairs of coupled two-state systems driven by a bichromatic external 
field is studied. We use a system of two coupled spin-i that can be translated into a three-state 
ladder model whose intermediate state represents the entangled state. We show that this entangled 
state can be prepared in a robust way with appropriate fields. Their frequencies and envelopes are 
■ derived from the topological properties of the model. 
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Entanglement is a key concept in various contemporary areas of active research in quantum physics. It explicitly 
demonstrates the non-local character of quantum theory, having potential applications in quantum communication, 
\Q \ cryptography and computation [Q . The preparation of an entangled state is of great interest for both fundamental 
£SJ ■ and applied reasons. During the last few years various methods for preparation of entangled states of atomic systems 
have been proposed and some of them experimentally demonstrated |^,|| . 

Although a quantum system can be manipulated by tailored sequences of resonant pulses of precise area, in particular 
ir and 7r/2 pulses respectively for the complete inversion and the equal weight coherent superposition, deviations from 
the precise pulse area and from resonance can lead to significant errors. Adiabatic passage techniques provide much 
■-"^ ■ greater robustness against fluctuations in the interaction parameters. The stimulated Raman adiabatic passage 
(STIRAP) method H has been proposed for the creation of an entangled state of two three-level atoms in a QED 
cavity || and for A atomic systems |J. 
■y ' In this paper we propose a simple method for entangling two subsystems driven by pulse shaped external fields. 
For definiteness we take these to be two identical spins interacting with each other and driven by radio-frequency 
fields. This system can be translated in a three-level ladder model with the intermediate level corresponding to the 
q-i entangled state 0. The goal is to populate completely this entangled state at the end of the pulses. The most efficient 
couplings are obtained with two near one-photon resonant fields. We show that unlike in the STIRAP process, one- 
and two-photon detunings are required to populate most efficiently the intermediate level. 



II. THE MODEL: TWO SPIN SYSTEM IN EXTERNAL FIELDS 



We consider two spin- ^ particles coupled by uniaxial exchange interaction in the z direction [[7J. In a time-dependent 
magnetic field B(f) = [B x (t), B y (t), B z ], the Hamiltonian of this system reads (ft = 1) 

H(t) = H + /xB(t). (§i + S 2 ), (1) 

where 

Ho = 4£Sf (8 Sf (2) 

is the part describing the exchange interaction, /i is the gyromagnetic ratio, £ is the exchange- interaction constant, and 
Sfc = [S%, S%, 5|] is the kth spin operator (k = 1,2). We assume that the static magnetic B z field in the z direction 
is strong enough (\nB z \ 3> |£|) so that the model with exchange interaction only in the z direction is justified. 

In the spin product state space {|m) 1 \m} 2 } (m =1,1), where the states \i) k and denote respectively the 
spin-down and spin-up states of the fcth spin, a complete basis of orthonormalized eigenstates of Hq is given by 
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Ill) = U>1 U> 2 , 

|IT + ) = ^ HOi IT> a + |T>! U> 3 ] , 

ITT) = 1T>! IT> 2 , 
|ir) ^[IDxIDs-lThlDJ. 



In this basis, the Hamiltonian (Q) can be expressed in the block-matrix form 

H(t) 

where 



H c (t) 



(3a) 
(3b) 
(3c) 
(3d) 



(4) 











(5) 



with /? = \fl x ,p y ,p z ] — /iB. The state ||T _ ) i s thus decoupled from the other states; it describes the evolution of 
a spin— singlet in a time-dependent magnetic field. This decoupling justifies our choice of the basis. The three 
other states |TT + ) an d ITT) are coupled by the transverse (xy) magnetic field. To complete the definition of the 
problem, we suppose that initially the two-spin system is in the unentangled state |||). Our goal is to establish the 
conditions leading to the most efficient robust transfer into the entangled state ||T + )- 

We consider the case when the spin system interacts with a constant magnetic field in the z direction and two 
radio-frequency fields of respective frequencies u\ and lo 2 in the x direction: 



P z = const, 

p x = ili (t) cos (ujit + ipi) + Q 2 (t) cos (to 2 t + fa) , 

Py=0. 



(6a) 
(6b) 
(6c) 



The state vector <fr(t) is solution of the Schrodinger equation i4f(p(t) — H c (t)</>(i) with the Hamiltonian H c (t) (j|) 
written in the basis {|||) , |IT + ) ; ITT)}- When the radiofrequency fields are off (P x = 0), we have the following 
energies En = £ — p z , E^+ = — £ and E^ = £ + p z . Without loss of generality we assume £ < and p z > 0, 
leading for p z > 2 |£| to a ladder configuration E±± < E^+ < E^. We apply near resonant fields u% « — En, 

u>2 ~ -Eft — ^1T+' i- e - w ith the detunings Ai and A2: 



u 1 = 2\£\+p z -A 1 , 
Lo 2 = -2\£\+p z -A 2 . 



(7a) 
(7b) 



According to the rotating wave approximation (RWA), one can neglect non resonant counterrotating terms under the 
conditions u>i, 2 3> S7i (t) , Q 2 (t). The RWA transformation 



e '" u 





L Q 










-i(E ll +uj 1 +w 2 )t 



leads to the state vector <p(t) = RJ(f>(t) (whose coefficients have the same absolute values as the ones of 
satisfies the Schrodinger equation 



(8) 



that 



i-</>(t) = H e (tMt), 

with the Hamiltonian H c = RtH c R — i(dR^ /dt) R, where only the quasi-resonant terms have been kept: 



(9) 



Oi (t) + e i5t n 2 (f) 



fil (t) 



-i6i 



n 2 (t) 



2Ai 



n 2 (t) 



,iSl 



fil (t) 



V 2 (t) + e-" 5 ^! (i) 2(Ax+A 2 ) 



(10) 
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The frequency 



6 = oji - u 2 = 4 |£| + A 2 - Ai 



(11) 



characterizes the coupling ambiguity ||. This Hamiltonian H c allows indeed both fields to couple the two transitions 
when the field amplitudes fli (t) and f2 2 (t) are not small compared to 5. The competing coupling schemes are 
depicted in Fig. pi with |1) = |jf) , |2) = ||f + ) an( i |3) = Iff). Two limit channels can thus be exhibited: the channel 
A (respectively B) corresponding to the absorptions of one u>\— photon (respectively of one w 2 — photon) by the 1 — 2 
transition and of one u> 2 — photon (respectively of one u)\— photon) by the 2 — 3 transition. We will study more precisely 
in Sections V and VI the various regimes that occur in this system. The problem of preparing the entangled state 
|2) = ||f + ) is t nus reduced to the study of the population transfer of the intermediate level in the ladder system 
driven by the Hamiltonian H c (|To|). 

The populations given by the Schrodinger equation (^|) are invariant under the following transformation T 



Ai -> Ai + 6, 
A 2 -> A 2 - 5, 
5^-5, 



(12a) 
(12b) 
(12c) 
(12d) 



We indeed obtain Rt (Yh c ) R-i(9Rt/dt) R = H c , with the unitary transformation 
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e 



-i6i 







(13) 
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III. NUMERICAL RESULTS 

Figures || and || display the population of the state |2) at the end of a sequence of delayed gaussian pulses of the 
same lengths and the same peak amplitudes 

fii (i) = Q a exp [- (t + r f /T 2 ] , (14a) 

Q 2 (i) = n exp -(t- t) 2 /T 2 , (14b) 

for various normalized peak amplitudes CIq/8 and detunings A/S, where we have chosen 

A = A 1 = A 2 . (15) 

We have considered this restriction of the parameters because it gives preferentially large islands of good population 
transfer. This will be justified in Sections V and VI. Note that the case Ai = — A 2 is irrelevant since it corresponds 
to a two-photon resonance between the product states \H) and Iff). We could have considered equivalently the 
restriction A 2 = Ai + 25 in accordance with the symmetry (fi2|). The two different orderings of pulses have been 
considered: the sequence 1 of Fig. || (respectively the sequence 2 of Fig. [|) corresponds to the wi— pulse (respectively 
the u>2— pulse) being switched on first, with the delay r = 1.7T (respectively r = — 1.7T). Global adiabaticity is 
ensured by the choice of a large pulse area ^IqT — 50. 

One can distinguish three islands of robust high transfer (white regions). Specific parameters characterizing these 
islands are labelled by (a), (a'), (d), and (df), with the subscript 1 or 2 respectively for the figures || or || [except (a 2 ) 
which is outside the regions of high transfer]. These islands of high transfer are analyzed in the following Sections 
by using the dressed Hamiltonian corresponding to H c ( |To| ) and the adiabatic properties of the dynamics. We will 
characterize different regimes and associate them with different effective dressed Hamiltonians. We will show that 
the islands of good transfer can be understood from the topological properties of the appropriate effective dressed 
Hamiltonian. 

We will show the following results: 

(i) Regions (a) correspond to a STIRAP-like process associated to the channel A (see Fig. |l|) that is perturbed (in 
the sense of non-resonant perturbation theory) by the channel B. Note that the restriction A 2 = Ai + 26 would have 
given a STIRAP-like process associated to the channel B perturbed by the channel A; 
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(ii) Regions (d) (in the weak field regime, i.e. Oi , < 8) correspond to an effective two-level SCRAP-like (Stark 
chirped rapid adiabatic passage) process JUnJ ; 

(iii) Regions (d') (in the strong field regime, i.e. fli,il 2 ^ $) corresponds to an effective two-level bichromatic 
SCRAP process (with additional Stark shifts) as described in Ref. [O. 



IV. THE DRESSED HAMILTONIAN 



It is convenient to use the adiabatic Floquet theory in order to study the Hamiltonian H c ( p^| ) since its time depen- 
dence contains a characteristic frequency S. The dressed Hamiltonian (or quasienergy Hamiltonian) corresponding to 
H r is SOl: 



K in u n 2] = _ i§ d + 1 
08 2 



fti (t) + e- w n 2 (t) 

Ox (t) + e w n 2 (t) 2Ai fl 2 (*) + e ie ^i (f) 

Q 2 (t) + e^Oi (t) 2 (Ai + A 2 ) 



(16) 



This dressed Hamiltonian depends parametrically on the pulse shapes and the detunin gs. It acts in the Hibert 
space spanned on the three states {|1) , |2) , |3)} tensored by the photonic Hilbert space The eigenstates of 

K are families of three states denoted |1; k, —k), |2; k — 1, —k) and |3; k — 1, — k — 1) with k a positive or negative 
integer. The corresponding eigenvalues Ai ; fc,-fc, M-,k-i,-k and \3- t k-i.-k-i have the following periodicity property: 
\i-kiM = ^n;fci— i,fca+i + ^1 f° r n — 1,2,3. The notation |n;fci, k 2 ) characterizes (when the fields are off) the state 
\n) dressed by the field of k\ lo±— photons and of k 2 lo 2 — photons. The integers k\ and k 2 characterize thus relative 
photon numbers of the respective fields of frequency u)\ and u 2 . The initial state is denoted |1; 0, 0) . The problem can 
be formulated as follows: we look for robust adiabatic connections between the initial state |1; 0, 0) and the final state 
|2; k — 1, —k) for some positive or negative integer k. 

The possible connections depend on the topology of the eigenenergy surfaces of ( fHf ) as functions of the field envelopes 
fii and n 2 for given Ai and A 2 . The topolo gy is characterized by true crossings which occur generically 

when one of the fields is off. We will study in the following the topology of K using different effective dressed 
Hamiltonians corresponding to different regimes. These regimes will depend on the ranges of the detunings and of 
the field amplitudes. 

We classify the different regimes and construct effective dressed Hamiltonians by determining in the Hamilonian K 
( |l6| ) which terms are resonant (or quasi-resonant) and which one are only perturbative. The resonant terms are treated 
by an adapted unitary transformation which allows an explicit diagonalization, whereas the perturbative terms can 
be treated by stationary pertubation theory. This technique has been presented in (lj] . Note that for a simple RWA 
two-level system of Rabi frequency f2 and detuning A, the perturbative regime is such that fl <C |A| and the resonant 
regime such that f2 > |A|. We classify the different regimes as functions of the ranges of the field amplitudes and of 
the detunings. In the following, we have normalized all the quantities with respect to 6. 



V. WEAK FIELD REGIME 



The weak field regime occurs when fii (t) , 2 (t) < 5. Note that when one has Ai = A 2 additionally, this regime 
coincides with a strong spin coupling since we have then 4 |£| > f^i (£) , fl 2 (t). In this case of weak field regime, we 
can intuitively analyze the different regimes with respect to the range of the detunings using the diagram of linkage 
patterns (Fig. |l|). Five relevant regimes (bounded by dashed lines) have been collected in Fig. |[ depending on the 
quasi-resonances : 

(A) The transition 1-2 is quasi-resonant with u>x and perturbed by oj 2 , 2-3 is quasi-resonant with lu 2 and perturbed 
by wi; 

(B) 1-2 is quasi-resonant with lo 2 and perturbed by uj\, 2-3 is quasi-resonant with lo± and perturbed by lo 2 ; 

(C) 1-2 and 2-3 are both quasi-resonant with uj 2 , and perturbed by u>x; 

(D) 1-2 is quasi- resonant with u> 2 and perturbed by lo\, 2-3 is perturbed by oj 2 and wi; 

(C) 1-2 and 2-3 are both quasi-resonant with u>x, and perturbed by lo 2 ; 

(D) 1-2 is quasi- resonant with u>\ and perturbed by lo 2 , 2-3 is perturbed by uj\ and u> 2 . 
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In the exact resonant cases, we have represented the regimes A, C and D in Fig. (g). 
As shown schematically in Fig. |[ the above regimes can be roughly bounded by 

Ai = ±5/2, Ai = -35/2 A 2 = ±6/2, A 2 = 35/2. (17) 

By the symmetry (|l2|), we recover the regime B from the regime A, C from C, _D from D (exchanging additionally 
fli and f2 2 ). 

The regimes A and B are STIRAP-like regimes; I? and D are SCRAP-like regimes. 

We do not consider other regimes where the state 1 1) is almost not depopulated by adiabatic passage. 

The line A = —6/2 appears as a dashed line in Figs g and ||, where the restriction Ai = A 2 has been considered. 



A. Regime A 



When the transition 1-2 is quasi-resonant with the frequency ui\ and the transition 2-3 quasi-resonant with the 
frequency lo 2 , the process can be analyzed as the channel A 'perturbed (in the sense non-resonant perturbation theory) 
by the channel B. We refer to it as the regime A as shown in Figs || and [| where it is roughly bounded by the dashed 
lines A = —5/2, A = 5/2 (not shown) and f^o = 5. This regime is approximately characterized by the following 
effective Hamiltonian in the basis {|1; 0, 0),|2; -1, 0),|3; -1, -1)} @: 



[»2(t)] 2 

"2(<5+Ai) 

fti (t) 2Ai 




fii (t) 

[^2(t)] 2 , 

2(«+Ai) 

n 2 (t) 



2(5-A 2 



2(A X 





n 2 (t) 

A 2 )- 



[Qi(t)] 2 

2(«-A 3 ) 



(18) 



which corresponds to the Hamiltonian characterizing the channel A with additional time dependent Stark shifts (on 
the diagonal) induced by the channel B. Note that this effective Hamiltonian is less precise for bigger or fl 2 
approaching 5. 

Before analyzing the dynamics given by this Hamiltonian ( |l8| ) , we recall the results in the limit case of a very weak 
field Sli (t) , f2 2 (t) <C 5 obtained in Ref. p3[ . In this case the perturbative terms can be neglected and the Hamiltonian 
becomes 



fii (t) 

n x (t) 2A1 n 2 (t) 

n 2 (t) 2(A! + A 2 ) 



(19) 



This resulting effective Hamiltonian corresponds to the channel A alone. The topology of the energy surfaces of this 
Hamiltonian has been analyzed in It has been shown that the adiabatic transfer to state |2) is topologically 

allowed for 



AiA 2 > 0. 



The topological analysis shows moreover that for the sequence 1 the region of this process is bounded in the parameter 
space by the curves 

fto = 2VAi(A! + A 2 ) 

and for the sequence 2 by the curves 

n Q = 2^Ai (Ai + A 2 ) and Q = 2^A 2 (Ai + A 2 ). 

Taking now into account the perturbation by the channel B [Hamiltonian (|lj)] leads to two kinds of topology as shown 
in Figs D an d @, where the surfaces of quasienergies as functions of the normalized Rabi frequencies £l\/5 and Q 2 /5 
respectively for A = Ai = A 2 = —5/20 and A = Ai = A 2 = —5/4 have been displayed. The eigenvalues of ( |l8| ) (not 
shown) fit these surfaces well except in Fig. |?j when f^i ~ 5 and il 2 ~ 5 because of an additional dynamical resonance 
(i.e. a resonance occuring beyond a threshold of the field amplitudes) [fl4"p|] which involves the surface connected with 
|3; 0, —2) (which corresponds to the surface connected to |3; —1, —1) and translated of 5) and the surface right below. 

Fig. ^ shows that the two conical intersections, one occuring for Qi = 0, the other one for fl 2 = 0, bound the 
adiabatic connection between the initial state |1; 0, 0) and the target state |2; —1,0). More precisely, for the sequence 
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1, the conical intersection occuring for f2j = is favorable for this adiabatic connectivity. The path denoted a\ (also 
corresponding to the cross a\ of Fig. |2|) is an example for the complete transfer. On the other hand, for the sequence 

2, the conical intersection occuring for fii = is also favorable but the one occuring for = is detrimental since it 
makes |1; 0, 0) connect with |3; — 1, —1). The path denoted a 2 is an example for the complete transfer to |3; — 1, —1) 
(also corresponding to the cross a 2 of Fig. ||). 

For a bigger detuning (in absolute value), the topology is different as shown in Fig. 0. The previous conical 
intersection occuring for Q 2 = has now disappeared and another one involving the surfaces connected to 1 1 ; 0, 0) 
and |3; 0, —2) has appeared. The two conical intersections, the one occuring for Six = and the other one for = 0, 
are involved for the adiabatic connection between the initial state |1;0, 0) and the target state 1 2; — 1,0). More 
precisely, for the sequence 1, these two conical intersections bound this adiabatic connection; the path denoted a[ 
(also corresponding to the cross a[ of Fig. ||) is an example for the complete transfer. On the other hand, for the 
sequence 2 only the conical intersection occuring for Oi = bounds now the adiabatic connection; the path denoted 
a' 2 is an example for the complete transfer (also corresponding to the cross a' 2 of Fig. ||). 

Using the effective Hamiltonian (|l^), the position of the previous conical intersections, for Qi = and O2 = 
respectively, lead to the three boundaries for the sequence 1 



A 



16(5 



-5^ 2 ± J 9 (a 



325 2 



Oi = W2(<5 - A) 2(8 + A) — y/2(S + A) 



(20a) 



(20b) 



The delay between the pulses has been chosen sufficiently large such that it is a good approximation to consider that 
the adiabatic connectivity is quite well described by the value of the peak ampl itude s. Thus we have displayed these 
boundaries in Fig. || as full lines, with f2 = ^ 2 for (20a) and with Sl = fii for ( |20b| ). They globally bound the lower 
and upper part of the island of good transfer of the regime A observed in the numerical computation. This island 
is crossed by the line of resonance A = around which the transfer to |2) depends on the pulse areas, as shown by 
small oscillating islands. 

For the sequence 2, the conical intersections involved give the following boundaries 



166 



9(n 2 f + 32S 2 



Qi = J(5-A) 4A + S ± y/S (6 + 8A) 



for A < 0. 



(21a) 



(21b) 



These curves are displayed in fig. g, with flo = fl 2 for (21a) and with Slo = ^1 for (21b). They give a good prediction 
of the island of good transfer of the regime A observed numerically. 

For the two sequences, the islands of good transfer to the state |2) of the regime A occur with absorption of one 
photon. 



B. Regime B 



This regime is characterized by the transition 1-2 quasi-resonant with the frequency 0J2 and the transition 2-3 quasi- 
resonant with the frequency w%. This process can be analyzed as the channel B perturbed (in the sense of non- resonant 
perturbation theory) by the channel A and is described by the effective Hamiltonian = TH^: 



2Ai 

n 2 (t) 




2(Ai +5) 



n 2 (t) 

[Q!(t)1 2 

r 2Ai 
fti (t) 



2A 2 



2 (Ai 





■ A a )+ [ ° 2(t)]2 



2A 2 



(22) 



The regions of high transfer efficiency to the state |2) are bounded in the same manner as in the regime A by the lines 
( pcf ) and (|2l]) to which we apply the transformation T (|l^) . 
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C. Regime C 



The regime C is characterized by a mixture of regimes A and B for which the transitions 1-2 and 2-3 are both 
quasi-resonant with the same frequency uj 2 . As long as the field is perturbative for both transitions, we have the 
following effective Hamiltonian 



H 



n 2 (t) 2 (Ai + s) + 







2Ai 

Q 2 (i) 



PM*)] 2 

2(A 2 -5) 



2 (Ai 





n 2 (t) 

A 2 + J) - 



PM*)] 2 



2(A 2 -5) 

in the basis {|1; 0, 0),|2; 0, —1) , |3; 0, —2)}. No efficient transfer is observed in this regime. 



(23) 



D. Regime D 



This regime is such that the only quasiresonance is between the states 1 and 2 with u> 2 . In this case, in the basis 
{|1; 0, 0) , 1 2 ; 0, — 1),|3; 0, —2)} we can construct an effective Hamiltonian from the previous one [Eq. (|23f)1 considering 
that the uj 2 — field is perturbative for the transistion 2-3: 



2 



2Ai 

n 2 (t) 2 (Ai + S) 




n 2 (i) 

[^(t)] 2 [Qi(t)] 2 



2Ai 2(A 2 -5) 




[»2(t)] 2 

2A 2 



2 (Ai 



Pl(t)] 2 

2(A 2 -«5) 



[»2(t)] 2 
2A 2 



(24) 



We can remark that this Hamiltonian is valid for the field amplitude f2 2 below the position of the resonance occuring 
between the transition 2-3 and the uj 2 — field that can be estimated by 



n T 2 = 2y/A 2 (Ai + A 2 + 5) and Ai + 2A 2 + S < 0. 



(25) 



This limit is represented as the bent dashed line crossing the figure horizontally in Figs || and || (with Qq — ^2) ■ Below 
this limit, one is allowed to decouple the states |2;0, — 1) and |3;0, — 2) from the Hamiltonian (|23|). A more detailed 
analysis of this regime shows that a dynamical resonance between the transition 1-2 and the tt>i— field, induced by the 
u) 2 — field occurs approximately for 



n 2 = nf = V-Ai (aT+25). 



(26) 



It is obtained when the difference of the dressed eigenvalues connected to |1;0, 0) and |2;0, — 1) (calculated without 
the Stark shifts) compensates the difference of the frequencies 5. This additional resonance is described as dynamical 
since it occurs beyond a threshold of the u> 2 — field amplitude. It is represented as the bent dashed line crossing the 
figure vertically (which separates the regimes D and D') in Figs || and || with Qo = il 2 r . The Hamiltonian ( p4| ) is 
thus approximately valid before the dynamical resonance (|2^) . 

Below this dynamical resonance, this Hamiltonian (24) is very similar to the one describing the process named 
SCRAP (acronym for Stark chirped Rapid adiabatic passage) between the states |1; 0, 0) and |2; 0, —1) ||. The pump 
of this process is here fi 2 and the Stark pulse f2i. We have here f2 2 acting additionally as a Stark pulse. 

It is important to note that when Ai = —8, the field uo 2 is exactly in resonance with the transition 1-2, and it cannot 
induce any complete population transfer from |1) to |2). Below this boundary (plotted as a full line in Figs ^ and 
||), i.e. for Ai < — 6, the topology does not allow the transfer from |1) to |2). Above this aboudary (Ai > —5), the 
transfer is possible as shown by the surfaces of quasienergies (for A = Ai = A 2 = 9(5/10) in Fig. ||. The eigenvalues 
of ([24j) (not shown) fit well these surfaces below the dynamical resonances fi 2 < f2 2 . Fig. || shows that the conical 
intersection for U 2 — between the surfaces connected to |1; 0, 0) and the target state |2; 0, —1) bound the adiabatic 
connection between these states. This characterizes a transfer to the state |2) with absorption of one w 2 — photon. 
This boundary is calculated from the effective Hamiltonian 



tlx = 2\ A 



S 2 



25 - A 



It is plotted in Figs |^ and || as a full line in the region D and bounds the upper island of good transfer of this region. 
The cases beyond the dynamical resonance are studied in the next section. 
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VI. STRONG FIELD REGIMES 



The strong field regime occurs when fii (t) , Q2 (t) > S. For Ai = A2, this corresponds to a weak spin coupling since 
one has then 4 |£| < f2i (t) , (*)• More resonances occur in this case and the previous effective Hamiitonians are no 
more vaiid. We wiil study in detail the interesting regime D' which gives quite large areas of transfer to state |2). 

This regime is located below the resonance (|5|) and beyond the dynamical resonance (^), when the transition 
1-2 is quasi- resonant with both the lo\— and ^—fields and when the transition 2-3 is not resonant with either the 
wi— field or the field. This regime is thus characterized by the effective dressed Hamiltonian 



-iS 



00 



n 2 (t) 











9 rji(t) 2(Ai 



n 2 (t) 

+ 6)- 



'i(> 



Slj. (t) 



2(A 2 -<5 





["2(t)] 2 

2A 2 



2 (A, 



6) 



2(A 2 -<5 



2A 2 



(27) 



It is equivalent to a two- level system driven by a bichromatic field |Tl[] with additional Stark shifts. The surfaces 
of quasienergies as functions of the normalized Rabi frequencies Sl±/6 and (for A = Ai = A2 = —78/5) are 

displayed in Fig. ^. This figure shows that the two conical intersections, one for fli — and one for CI2 = 0, bound 
the adiabatic connection between the initial state |1; 0, 0) and the target state |2; 1, — 2). We calculate the boundaries 
using the effective Hamiltonian (^7j) , which are plotted as full lines in Figs || and ||: 



2W(A-<5) 26- A -2^5 {25- A) 



2 = 2, A 



d-A-y/6(6 — 4A) 



This process corresponds to a multiphoton transfer to the state |2), with absorption of two loi— photons and emission 
of one u)\— photon. 

The analysis of the topology allows to improve the transfer efficiency. It shows indeed that a wi— field amplitude 
larger than the loi —field amplitude is better in this regime since the conical intersection for f^i — occurs for a larger 
value than the one for fi 2 = 0. 

This process of a two-level system driven by a bichromatic field studied in Ref |llj shows that the transfer can 
still occur for a stronger field (i.e. for a weaker spin coupling), but with absorption of more than two loi— photons 
and emission of more than one u)\— photon. This result is shown in Fig. [l^ where strong field white islands can be 
observed. The lower white islands correspond to good population transfer to the entangled state |2; k — 1, —k), with 
k = 1,2,3,4 from left to right. 



VII. CONCLUSION 



We have shown that a system of two interacting spins in an external bichromatic field is equivalent to a three 
level problem. We have determined the choices of laser pulses which can give a maximal final population in the 
entangled state. The proposed strategies are robust with respect to the external parameters. We have found that in 
the parameter space it is possible to find large islands where the quantum system can be transferred to the entangled 
state with a high efficiency. These islands of good transfer have been characterized by the topology of the surfaces of 
dressed states as functions of the parameters. 

The methods employed here are quite general and can be applied for a large variety of systems. We anticipate 
interesting applications of this method in quantum computing and quantum communication. 
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FIG. 1. Diagram of linkage patterns between the three states showing the different couplings. Note that Ai and A2 have 
been chosen here negative. 
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FIG. 2. Contour map of population transfer efficiency P2(oo) for varying peak Rabi frequency fio/5 and varying detuning 
A/5 (with A = Ai = A2) for the sequence 1. White areas correspond to high efficiency transfer (close to 1) to the entangled 
state. Dark areas correspond to low efficiency transfer (close to 0) to the entangled state |2). The dashed lines separate different 
regions labelled A, D and D', associated to different effective Hamiltonians constructed in Sections |y] and |vj|. The regimes of 
good population transfer are bounded by full lines predicted from the topological analysis. The crosses labelled (ai), (ai), (di), 
and (di) refer to parameters leading to high efficiency. They also refer to the pathways shown respectively in Figs ^, || and 
^| The regions A and D 

FIG. 3. Contour map of population transfer efficiency P2(oo) for varying peak Rabi frequency fio/5 and varying detuning 
A/5 for the sequence 2. The cross labelled (a' 2 ) in one region of low efficiency and the ones labelled (02), (C2) and (efo) in the 
regions of high efficiency refer to the pathways shown respectively in Figs (?], || ^| and 9. 

FIG. 4. Schematic diagram of the regimes for a weak field regime as a function of the normalized detunings A1/8 and A2<5. 
The restriction A2 = Ai has been used for Figs. H and H. 



FIG. 5. Diagram of linkage patterns for the three regimes (in the resonant case): A (Ai = A2 = 0), C (Ai = —5, A2 = 0) 
and D (A a = -5, A 2 = -5). 

FIG. 6. Quasienergy surfaces (in units of 5) as functions of fii/5 and Q2/8 for Ai = A2 = —5/20. The path denoted a\ 
(sequence 1), for fio = 0.355, connect the states |1) and |2) with the absorption of one oj\— photon. The path denoted 02 
(sequence 2), for fio — 0.355 connect the states |1) and |3) with the absorptions of one uii— photon and of one u>2 — photon. 

FIG. 7. Quasienergy surfaces as functions of fii/5 and Q2/S for Ai = A2 = —8/4. The two different paths, denoted a\ and 
02 (for fio = 0.75) depending on the sequence of the pulses connect the states |1) and |2) with the absorption of one wi- photon. 

FIG. 8. Quasienergy surfaces as functions of fii/5 and Q2/S for Ai = A2 = —95/10. Two different paths (denoted d\ and 
0I2) for fio = 0.85 connect the states |1) and |2) with the absorption of one UJ2— photon. 

FIG. 9. Quasienergy surfaces as functions of Qi/5 and Q.2/8 for Ai = A2 = —75/5. Two different paths (denoted d' x and 
d' 2 ) for fio = 35/2 connect the states |1) and |2) with the absorption of two u>2— photon and the emission of one wi— photon. 



FIG. 10. Contour map of population transfer efficiency P2(oo) as in Fig. @, but for stronger field amplitudes. 
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